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(~| , Combined data from gamma-ray telescopes and cosmic-ray detectors have produced 

^ J some new surprising insights regarding intergalactic and galactic magnetic fields, as well 

as extragalactic background light. We review some recent advances, including a theory 
explaining the hard spectra of distant blazars and the measurements of intergalactic 
magnetic fields based on the spectra of distant sources. Furthermore, we discuss the 
possible contribution of transient galactic sources, such as past gamma-ray bursts and 
hypernova explosions in the Milky Way, to the observed flux of ultrahigh-energy cosmic- 
rays nuclei. The need for a holistic treatment of gamma rays, cosmic rays, and magnetic 
fields serves as a unifying theme for these seemingly unrelated phenomena. 
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1. Gamma ray astronomy of cosmic rays 

Gamma rays from Active Galactic Nuclei (AGN) are studied extensively using 
ground-based atmospheric Cherenkov telescopes (ACT), as well as Fermi Space 
Telescope and other instruments. Their signals reveal important information about 
the sources, as well as about extragalactic background light (EBL) and intergalactic 
magnetic fields (IGMF) along the line of sight. The same sources are expected to 
accelerate cosmic rays, although it is more difficult to associate cosmic rays with 
their sources because the local, galactic magnetic fields alter the arrival directions 
of cosmic rays. 



1.1. Secondary gamma rays from the line-of -sight interactions of 
cosmic rays 

It was recently proposed that the hardness (and uniform redshift-dependent shape) 
of gamma-ray spectra of distant blazars can be naturally explained by the line- 
of-sight interactions of cosmic rays accelerated in the blazar jets.'^^ The cosmic 
rays with energies below 10^^ — 10^^ eV can cross large distances with little loss of 
energy and can generate high-energy gamma rays in their interactions with cosmic 
background photons relatively close to the observer. Such secondary gamma rays 
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can reach the observer even if their energies are well above TeV. In the absence 
of cosmic-ray contribution, some unusually hard intrinsic spectrgP or hypothetical 
new particleJiSlIil have been invoked to explain the data. 

As long as the IGMFs are smaller than '^10 femtogauss, secondary gamma rays 
come to dominate the signal from a sufficiently distant source. One can see this 
from the way the flux scales with distance for primary and secondary gamma raysP 
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Obviously, for a sufficiently distant source, secondary gamma rays must domi- 
nate because they don't suffer from the exponential suppression as in Eq. ((TJ. The 
predicted spectrum turns out to be similar for all the distant AGN, depending only 
on their redshift. These predictions are in excellent agreement with the data.l^tSl 

One can see the transi- 
tion from primary to secondary 
gamma rays in Fig. [TJ which 
shows the spectral index dif- 
ference for blazar spectra as a 
function of their redshifts. At 
small redshifts, the data con- 
firm the Stecker - Scully re- 
lation,'i2l but, at redshifts 0.15 
and beyond, there is clearly 
a new population of blazars, 
whose observed spectral in- 
dex shows only a weak de- 
pendence on the redshift. The 
nearby population is obviously 
the blazars from which primary 
gamma rays are observed. The 
distant blazars are observed in 
secondary gamma rays, which 
are produced in line-of-sight 
cosmic ray interactions. These 

secondary gamma rays are produced relatively close to the observer, regardless of 
the distance to the source. Hence, their redshift dependence is much weaker.!^ Fi- 
nally, there is an intermediate population around redshift 1.2 which is composed of 
some blazars seen in primary gamma rays and some seen in secondary gamma rays. 




Fig. 1. Spectral index change <5r = FqcV — FtcV as a 
function of redshift. While the low-redshift blazars agree 
with the Stecker— Scully relation,C21 the data indicate the 
existence of an additional, distinct population with a weak 
redshift dependence at redshifts 0.15 and beyond."" In par- 
ticular, the recently measured redshift^^ of PKS 0447-439 
is in agreement with the trend. 
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A recent redshift measurement of PKS 0447-439 redshiflP^ further strengthens 
our interpretation. Gamma rays with energies above 1 TeV have been observed 
from this blazar by HESS.^** The spectral properties agree with the trend (Fig. [Ij. 
Furthermore, there is no way for primary gamma rays to reach Earth from such a 
distant source, while secondary gamma rays provide a consistent explanation of the 
PKS 0447-439 spectrum.E^ 

This motivates future observations by ACT of blazars with known large redshifts. 
Secondary gamma rays with TeV and higher energies can be observed even from 
some sources located at cosmological {z ^ 1) distances.!^ 

The spectral slope of pro- 
tons and the level of EBL do 
not have a strong effect on 
the spectrum of secondary pho- 
tons, as one can see from Fig. [2] 
However, for the same photon 
flux, the neutrino flux varies 
depending on the maximal en- 
ergy -©max to which the protons 
are accelerated. Indeed, there 
are two competing processes 
that generate secondary pho- 
tons: pjEBL V^^ ^ Pll and 
PlCMB — >■ pe'^e~ . For smaller 
-E'max, a larger fraction of pho- 
tons come from the hadronic 
channel, which is accompanied 
by production of neutrinos via 

PJEBL JiTT* followed by the decays of charged pious and the neutron. Neutrino 
observations can help determine this parameter.^^l 



1.2. IGMFs and EBL 

The success of this picture lends support to the hypothesis of cosmic ray acceleration 
in AGN. Furthermore, one can use the spectral gamma-ray data to study EBL and 
IGMFs. The predicted spectra depend to some extent on the EBL model, as shown 
in Fig. [31 although this dependence is too weak to distinguish between different 
models.'^ IGMFs, however, can have a strong effect on the goodness of fit. Based on 
the spectra of several distant blazars, one can set both upper and lower limits on 

igmfS 
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Fig. 2. Photon (low energy) and neutrino (higli en- 
ergy) spectri2i expected from an AGN at ^ = 0.14 
(such as 1ES0229-I-200), normalized to HESS data points 
(shown),Ell for Emax = lO^GeV, 10i°GeV, and 10"GeV 
shown by the solid, dashed, and dash-dotted lines, respec- 
tively. 



10-i^G < B < 3 X 10-^^G. 
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Fig. 3. Sensitivity of secondary gamma-ray spectra to the model of EBL and to the average 
value of IGMFs. Fermi upper limit and tentative detectiorli^ for blazar lES 0229+200 are shown 
at lower energy, and HESS data points at high energy. The predictions of two models of EBL are 
shown in the left panel according to Essey et al. Eland the effects of intergalactic magnetic fields 
are shown in the right panel.El 



1.3. Time variability 



10 
10 
10 



14 



S 10^3 

5-10^^ 
S 10^° 
i 10^ 



An important property of secondary gamma rays is the lack on short-scale time 
variability.* For E > I TeV and z > 0.15, one expects the signal to be dominated 
by secondary photons, and any time variability on short scales should be erased by 
delays in the propagation of protons and electromagnetic cascades. Fig. |4] shows the 
time delays as a function of the proton energy. 

The present data are not 
yet sufhcient to probe time 
variability at the relevant ener- 
gies and redshifts because the 
data points are too few. While 
time variability has been ob- 
served for nearby TeV blazars 
at TeV energies, as well as for 
distant TeV blazars at energies 
above a few hundred GeV, no 
variability has been reported 
so far for distant TeV blazars 
at TeV energies. One can in- 
fer from Fig. [Dhow distant the 
source has to be for the sec- 
ondary signals to dominate. It 
is evident that the secondary 
component takes over for red- 
shifts beyond 0.15. 
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4. Time delays of gamma rays emitted at redshift 
2 = 0.17 for different proton energies Eq in a femtogauss 
random IGMF with a correlation length of 1 Mpc. 
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2. Composition of UHECR and past transient phenomena in the 
Milky Way 

Let us now turn to another phenomenon related to cosmic rays and magnetic fields, 
only this time we will concentrate on the magnetic fields inside the galaxy and their 
effect on the observed fluxes of ultrahigh-energy nuclei.GsHH There IS a growing 
evidence that long GRBs are caused by a relatively rare type(s) of supernovae, 
while the short GRBs probably result from the coalescence of neutron stars with 
neutron stars or black holes. Compact star mergers undoubtedly take place in the 
Milky Way, and therefore short GRBs should occur in our Galaxy. 

Although there is some correlation of long GRBs with star-forming metal-poor 
galaxies, many long GRBs are observed in high-metallicity galaxies as wellj'^iHzl] 
and therefore one expects that long GRBs should occur in the Milky Way. Less pow- 
erful hypernovae, too weak to produce a GRB, but can still accelerate UHECR,'^^! 
with a substantial fraction of nuclei 

If the observed cosmic rays originate from past explosions in our own Galaxy, 
PAO results have a straightforward explanation.^^ 

GRBs have been proposed as the sources of extragalactic UHECR,'2SH3ll and they 
have also been considered as possible Galactic sources .I^SHII] It is believed that GRBs 
happen in the Milky Way at the rate of one per ^grb ~ 10'* - 10^ yeavs^^^ Such 
events have been linked to the observations of positrons .1^^^^ 

If local sources, such as past GRBs, hypernovae, and other stellar explosions in 
the Milky Way, produce a small fraction of heavy nuclei,!^ the observed fraction of 
UHE nuclei is greatly amplified by diffusion. This is because the galactic magnetic 
fields are strong enough to trap and contain nuclei but not protons with energies 
above EeV. This observation leads to a simple explanation of the composition trend 
observed by PAO. 

As illustrated in Fig. [5l diffusion depends on rigidity, and, therefore, the ob- 
served composition can be altered by diffusion.^® Changes in composition due 
to a magnetic fields have been discussed in connection with the spectral "knee" ^ 
and also for a transient source of UHECR. ^ The "knee" in the spectrum occurs 
at lower energies than those relevant PAO, and at higher energies the cosmic rays 
effectively probe the spectrum of magnetic fields on greater spatial scales, of the 
order of 0.1 kpc.l^ 

One can use a simple modeP^ to show how diffusion affects the observed spec- 
trum of the species "i" with different rigidities. Let us suppose that all species are 
produced with the same spectrum nf''"^^ — n^^'^'^^ oc E^"' at the source located in the 
center of the Milky Way and examine the observed spectra altered by the energy 
dependent diffusion and by the trapping in the Galactic fields. 

In diffusive approximation, the transport inside the Galaxy can be described by 
the equation: 
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Fig. 5. For each species, there is a critical energy i?o,i for which the Larmor radius Ri is equal 
to the magnetic coherence length Ic- For E <^ Eg^i, the mean free path of the diffusing particle is 
I ^ lo, and Di{E) = lc/3. For E 3> Eq^, the particle is deflected only by a small angle 6 ^ lo/Ri, 
and, after k deflections, the mean deflection angle squared is ^ k{lQ/Ri)^. The corresponding 
diffusion coefficient is Di(E) oc (^gr^)^, for E 2> -Bo,i- 



Here Di{E, r, t) = Di{E) is the diffusion coefhcient, which we wiU assume to be 
constant in space and time. The energy losses and all the interactions that change 
the particle energies are given by hi{E) and the kernel in the collision integral 
Pik{E, E'). For energies below GZK cutoff, one can neglect the energy losses on the 
diffusion time scales. 

The diffusion coefficient D{E) depends primarily on the structure of the mag- 
netic fields in the Galaxy. Let us assume that the magnetic field structure is com- 
prised of uniform randomly oriented domains of radius Iq with a constant field B in 
each domain. The density of such domains is iV ~ Iq^. The Larmor radius depends 
on the particle energy E and its electric charge qi = eZi. 



The spatial energy spectrum of random magnetic fields inferred from observa- 
tions suggests that B ^ 3/iG on the 0.3 kpc spatial scales, and that there is a 
significant change at / = 1/fc ~ 0.1 — 0.5 kpc.^"* This can be understood theo- 
retically because the turbulent energy is injected into the interstellar medium by 
supernova explosions on the scales of order 0.1 kpc. This energy is transferred to 
smaller scales by direct cascade, and to larger scales by inverse cascade of magnetic 
helicity. Single-cell-size models favor '^0.1 kpc scales as welkin 





(4) 
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As explained in the caption of Fig. 1, diffusion occurs in two different regimes 
depending on wfietlier tlie Larmor radius is small or large in comparison with the 
correlation length. As a result, the diffusion coefficient changes its behavior dramat- 
ically at E = Eo^i'. 



D^{E) = 




E < Eo^i, 
E > Eo^i- 



(6) 



Here the two parameters < 61^2 < 0.5 are different from zero if the magnetic 
domains are not of the same size. The exact values of these parameters depend on 
the power spectrum of turbulent magnetic fields. 

The approximate solution of the transport equation in our simple model yields 
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Since diffusion depends on 
rigidity, the composition be- 
comes energy dependent. In- 
deed, at critical energy i?o,j, 
which is different for each nu- 
cleus, the solution ([7]) changes 
from oc E~'^ to oc E~'^~'^ be- 
cause of the change in Di{E), 
as discussed in the caption of 
Fig. 1. Since the change occurs 
at a rigidity-dependent critical 
energy £^o,i = eEoZi, the larger 
nuclei lag behind the lighter 
nuclei in terms of the critical 
energy and the change in slope. 
If protons dominate for E < 
Eq, their flux drops dramati- 
cally for E > Eq, and the heavier nuclei dominate the flux. The higher Zi, the 
higher is the energy at which the species experiences a drop in flux. 

The modeP^ provides a qualitative description of the data (see Fig.]!]). To repro- 
duce the data more accurately, it must be improved. First, one should use a more 
realistic source population model. Second, one should include the coherent compo- 
nent of the Galactic magnetic field. Third, one should not assume that UHECR 
comprise only two types of particles, and one should include a realistic distribution 
of nuclei. Finally, one should include the extragalactic component of UHECR pro- 
duced by distant sources, such as active galactic nuclei (AGN) and GRBs (outside 
the Milky Way). A recent realization that very high energy gamma rays observed 
by Cherenkov telescopes from distant blazars are likely to be secondary photons 
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UHECR spectra,E51for the magnetic field ~ lO^iG, 
coherent over Iq = 100 pc domains. The power and the iron 
fraction were adjusted to fit the PAO dataPSl 
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produced in cosmic ray interactions along the line of sight lends further support to 
the assumption that cosmic rays are copiously produced in AGN jetsP^For energies 
E > 3x 10^^ eV, the energy losses due to photodisintegration, pion production, pair 
production and interactions with interstellar medium become important and must 
be included. The propagation distance in the Galaxy exceeds 10 Mpc, so that the 
Galactic component should exhibit an analog of GZK suppression in the spectrum. 
Extragalactic propagation can also affect the composition around 10^^ eV.I^ 




Fig. 7. Monte Carlo simulation of a typical set of "hot spots" in the directions of the closest, 
recent GRBs or hypernovae, in the model of Calvez et a im The distribution of transient sources 
was modeled for long GRBs (left) and short GRBs (right). 

Galactocentric anisotropy for a source distribution that traces the stellar counts 
in the Milky Way is small.^^ Although the anisotropy in protons is large at high 
energies, their contribution to the total flux is small, so the total anisotropy was 
found to be < 10%, consistent with the observations. The latest GRBs do not 
introduce a large degree of anisotropy, as it would be in the case of UHE protons, 
but they can create "hot spots" and clusters of events (Fig. [T]). 

The modefiS] leads to the following prediction for the highest-energy cosmic rays. 
Just as the protons of the highest energies escape from our Galaxy, they should 
escape from the host galaxies of remote sources, such as AGN. Therefore, UHECR 
with E > 3 X 10^^ eV should correlate with the extragalactic sources. Moreover, 
these UHECR should be protons, not heavy nuclei, since the nuclei are trapped in 
the host galaxies. If and when the data will allow one to determine composition on a 
case-by-case basis, one can separate E > 3 x 10^^ eV events into protons and nuclei 
and observe that the protons correlate with the nearby AGN. This prediction is 
one of the non-trivial tests of our model: at the highest energies the proton fraction 
should exist and should correlate with known astrophysical sources, such as AGN. 
The microgauss magnetic fields in the Milky Way cause relatively small deflections 
for the highest-energy protons. As for the intergalactic magnetic fields, there are 
reasons to believe that they are relatively weak, of the order of a femtogauss,!^ 
and, therefore, they should not affect the protons significantly on their trajectories 
outside the clusters of galaxies. 

If local. Galactic GRBs are the sources of UHECRs, the energy output in cosmic 
rays should be of the order of 10^^ erg per GRB. This is a much lower value than 
what would be required of extragalactic GRBs to produce the same observable flux. 
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Indeed, in our model the local halo has a much higher density of UHECR than 
intergalactic space, and so the overall power per volume is much smaller. The much 
higher energy output required from cxtragalactic GRBJsMHI jn UHECR has been a 
long-standing problem. The same issue does not arise in our case because it seems 
quite reasonable that a hypernova or some other unusual supernova explosion would 
generate 10''*^ erg of UHECR with energies above 10 EeV. 



3. A gamma- ray signature of cosmic-ray nuclei 

A spectral feature, namely an "iron shoulder" at 5-10 GeV can help identify 
cosmic nuclear accelerators.!^ Nuclei are likely to come out of acceleration re- 
gions unstable because they can lose a nucleon or a few nucleons to photodisin- 
tegration in the high-density photon environments accompanying some accelera- 
tors.!^ An unstable nucleus decays, and most of such decays are /3-decays. With 
a probability of order one, the /3-decay electron is captured by the Coulomb 
potential of the fully ionized atom.1251 Hence, a non-negligible fraction of nu- 
clei come out of astrophysical accelerators in the form of one-electron ions. 
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In a narrow en- 
ergy range, CMB pho- 
tons have energies ~ 
7 keV in the rest frame 
of the ion. Such pho- 
tons can excite the ion, 
which later emits a 
7 keV photon (in the 
ion's rest frame). Mul- 
tiple excitations and 
de-excitations can take 
place resulting in emis- 
sion of gamma rays, 
which have energies 
of 5 - 10 GeV in 
the laboratory frame. 
The spectral feature 
around 8 GeV (Fig. ^ 
can be used for iden- 
tifying astrophysical sources 

of nuclei, or (in the case of non-detection) for setting the upper limits on nuclear 
acceleration.!^ 
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Fig. 8. Expected signature of nuclear emission for Cen A (solid 
line), normalized t o tota l observed flux of iron, and the data from 
FermiMland HESS.1^21511 
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4. Conclusions 

Based on the recent data, one can make several remarkable inferences about the 
ultrahigh-energy cosmic rays and magnetic fields inside the Milky Way and in the 
intergalactic space. Gamma-rays detected from most distant blazars are most likely 
dominated by the secondary photons produced in line-of-sight interactions of cos- 
mic rays. This interpretation allows one to set both upper and lower bounds on 
intergalactic magnetic fields, lO^^^G < S < 3 x IQ-^^gW 

P\irthermore, the energy dependent composition of UHECR, with heavier nuclei 
at high energy, points to a non-negligible contribution from Galactic sources.'^ 
Diffusion in turbulent Galactic magnetic field traps the nuclei more efficiently than 
protons, leading to an increase in the nuclear fraction up to the energy at which iron 
escapes (~ 30 EeV). At higher energies, the extragalactic protons should dominate 
the flux of UHECR, and their arrival directions should correlate with locations of 
the known sources. 

If and when the neutrino telescopes, such as IceCube,'^ detect point sources, 
one can learn about the cosmic-ray sources and photon backgrounds by compar- 
ing the neutrino flux to the photon flux. Neutrino and gamma-ray observations 
can help distinguish the local Galactic sources from extragalactic sources of UHE 
nuclei.'^^lt^ These inferences open exciting new opportunities for multi-messenger 
photon, charged-particle, and neutrino astronomy. 

This work was supported by DOE Grant DE-FG03-91ER40662. 
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